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Debido a cambios en el clima, en el uso del suelo, 
principalmente por urbanización, y a los aprovechamientos 
hidráulicos, los registros de datos hidrológicos extremos 
han dejado de ser estacionarios, es decir, sus propiedades 
estadísticas están cambiando en el tiempo. Cuando un 
registro de crecientes o lluvias máximas, ambas anuales, 
muestra tendencia estadísticamente diferente de cero, 
el análisis de frecuencias orientado a estimar sus eventos 
de diseño debe tomar en cuenta tal comportamiento no 
aleatorio. El modelo probabilístico más simple que se 
puede establecer para procesar registros no estacionarios 
emplea la distribución log-normal de dos parámetros de 
ajuste, con su media variando en relación con el tiempo 
como covariable explicativa, según una regresión lineal 
logarítmica. Este modelo permite estimar el impacto, por 
ejemplo, al final de una década a futuro, de la tendencia 
ascendente (descendente) en las predicciones obtenidas al 
final del registro, a través de dos indicadores: los factores de 
aumento (disminución) de las predicciones y de la reducción 
(aumento) de la recurrencia de un evento de un determinado 
periodo de retorno. El análisis de frecuencias no estacionario 
con tal modelo se aplica a dos registros de crecientes 
anuales de la Región Hidrológica núm. 10, Sinaloa, y uno de 
precipitación máxima diaria anual de la Región Hidrológica 
núm. 37, El Salado. Las conclusiones destacan la sencillez y 
utilidad del modelo descrito para abordar de forma inicial 
los análisis de frecuencias no estacionarios. 


Palabras clave: crecientes anuales, precipitación máxima 
diaria anual, distribución log-normal, regresión lineal, error 
estándar de ajuste, factores de aumento de las predicciones. 


Abstract 


Campos-Aranda, D. E. (May-June, 2016). Simple Probability 
Model for Frequencies Analysis in Extreme Hydrological Records 
with Trend. Water Technology and Sciences (in Spanish), 7(3), 
171-186. 


Due to climate changes, land use mainly by urbanization and 
hydraulic works, the records of extreme hydrological data are no 
longer stationary, that is that their statistical properties are changing 
over time. When an annual record of floods or maximum rainfall 
shows a trend statistically different from zero, the frequencies analysis 
aimed at estimating their design events should take into account 
such non-random behavior. The simplest probabilistic model that 
can be established to process non-stationary records relies on the log- 
normal distribution with two fitting parameters, with time-varying 
mean as explanatory covariate, according to a logarithmic linear 
regression. This model allows estimating the impact, for example at 
the end of a decade in the future, of the upward (dorwniward) trend 
on the obtained predictions at the end of a record. The latter through 
two indicators: the predictions magnification (diminution) factors 
and the reduction (increase) in the recurrence of an event of a given 
return period. Non-stationary frequencies analysis with such a model 
is applied to two records of annual floods of the Hydrological Region 
núm. 10 Sinaloa and one annual daily maximum precipitation of the 
Hydrological Region núm. 37 El Salado. Conclusions highlight the 
simplicity and utility of the described model as an initial approach to 
the non-stationary frequencies analyses. 


Keywords: Annual floods, annual maximum daily precipitation, 
log-normal distribution, linear regression, standard error of fit, 
predictions magnification factors. 
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Introducción y objetivo 


Cuando se dispone de un registro de datos 
hidrológicos extremos, como gastos, niveles, 
lluvias de cierta duración y vientos, todos ellos 
máximos anuales, se puede realizar un análisis 
de frecuencias para estimar predicciones de tales 
variables aleatorias, asociadas con determinadas 
probabilidades de excedencia, cuyo recíproco 
es el llamado periodo de retorno o intervalo 
promedio en años entre la ocurrencia de un 
evento igual o mayor. En resumen, el análisis 
de frecuencias permite estimar eventos de diseño 
y consta de cuatro pasos: (1) selección de un 
modelo probabilístico o función de distribu- 
ción de probabilidades acumuladas (FDP); (2) 
aplicación de uno o varios métodos de estima- 
ción de los parámetros de ajuste de la FDP; (3) 
adopción de la FDP que mejor representa al 
registro disponible, lo anterior, de acuerdo con 
un criterio de calidad estadística del ajuste, y 
(4) con base en la FDP seleccionada, estimación 
de las predicciones buscadas (Khaliq, Ouarda, 
Ondo, Gachon, éz Bobée, 2006). 

Cuando no existe el registro de crecientes 
requerido en el sitio de interés o es inadecua- 
do por amplitud escasa, se recurre al análisis 
regional de frecuencias, el cual se desarrolla en 
dos etapas: a) se busca y selecciona un grupo de 
cuencas aforadas, relativamente cercanas, con 
un régimen hidrológico similar al de la cuenca 
bajo estudio; b) se estiman las predicciones bus- 
cadas en el sitio sin datos, a través de diversos 
métodos de transferencia de información hidro- 
lógica de los sitios o registros identificados en 
la etapa A (Leclerc € Ouarda, 2007). La técnica 
también es aplicable a estaciones pluviográficas 
o pluviométricas. 

La validez de los resultados del análisis de 
frecuencias local o regional está basado, de ma- 
nera especial, en dos suposiciones clásicas del 
registro: que sus datos sean independientes y que 
sus propiedades estadísticas no cambien en el 
tiempo o condición de estacionario. En la práctica, 
la independencia se pierde cuando el muestreo 
de los datos es próximo o cercano, y entonces 
un valor puede estar ligado con el anterior. En 
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registros anuales, tal problema no ocurre. Lo no 
estacionario se origina por el cambio climático 
o por los efectos de las actividades humanas en 
la cuenca, principalmente los aprovechamientos 
hidráulicos, o bien los cambios de uso del suelo, 
sobre todo la urbanización (Khaliq et al., 2006; 
Vogel, Yaindl, € Walter, 2011). 

Por lo general, y con base en la zona geo- 
gráfica, el cambio climático puede generar 
tendencias sutiles ascendentes o descendentes 
en los registros de datos hidrológicos extremos. 
Sin embargo, las tendencias fuertes ascendentes 
observadas en los registros de crecientes anuales 
y de precipitación máxima diaria anual se origi- 
nan por la urbanización ocurrida en la cuenca 
o en las inmediaciones de la estación pluvio- 
métrica, y la tendencia severa descendente se 
asocia en los registros de crecientes con los 
aprovechamientos hidráulicos, y en las series 
de lluvias máximas con las deforestaciones o 
desecación de lagos, ambas cercanas (Vogel et 
al., 2011). 

El objetivo de este trabajo consistió en expo- 
ner con detalle el método probabilístico más 
simple que se puede establecer para el análisis 
de frecuencias en registros de datos hidrológicos 
extremos que presentan tendencia importante. 
Este modelo utiliza la FDP log-normal de dos 
parámetros de ajuste, en la cual su media varía 
linealmente con el tiempo (Vogel et al., 2011). Se 
aplica tal modelo al análisis de frecuencias de 
dos registros de crecientes anuales de la Región 
Hidrológica núm. 10, Sinaloa; el de la estación 
Badiraguato con tendencia ascendente y el de 
la estación Puente Subpacífico con tendencia 
descendente. "Iambién se realiza el análisis de 
frecuencias del registro de precipitación máxima 
diaria anual de la estación climatológica Zacate- 
cas de la Región Hidrológica núm. 37, El Salado, 
que presenta tendencia ascendente. 


Descripción del modelo probabilístico 
no estacionario 


FDP adoptada 


Como el objetivo del estudio consiste en iden- 
tificar la tendencia en los registros de datos 
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hidrológicos extremos anuales y estimar sus 
consecuencias en el análisis de frecuencias, es 
necesario utilizar una función de distribución de 
probabilidades (FDP) que describa la relación 
entre magnitudes de los datos y predicciones 
asociadas con determinadas probabilidades de 
no excedencia. En este caso, y de acuerdo con 
Vogel et al. (2011), se adopta la distribución log- 
normal de dos parámetros (LN2), que fue uno 
de los primeros modelos probabilísticos utiliza- 
dos en los análisis de frecuencias de crecientes 
y que aún sigue vigente con su versión de dos 
(Ashkar éz Aucoin, 2012) y tres parámetros de 
ajuste (Hosking éz Wallis, 1997). Sin embargo, la 
razón principal de tal selección radica en que la 
función LN2, cuando se combina con un modelo 
de tendencia logarítmico lineal, origina una FDP 
no estacionaria bastante simple, la cual permite 
estimar el impacto de la tendencia en la magni- 
tud y el periodo de retorno de las predicciones, 
en un lapso a futuro seleccionado. 

Un registro de datos hidrológicos extremos 
anuales es una serie anual de máximos x,, en donde 
t varía de 1 an su número de datos. Si tal serie es 
estacionaria, puede seguir una distribución LN2, 
cuya función de cuantiles que permite estimar 
las predicciones con probabilidad de excedencia 
p o de periodos de retorno T = 1/p es (Vogel et 
al., 2011): 


X, =exp(u, +2, :0,) (1) 


en la cual, u, y o, son la media y la desviación 
estándar de los logaritmos naturales de x, es 
decir, que y, = 1n(x,) y z, es el valor de la variable 
aleatoria normal estándar con probabilidad de 
excedencia p. 


FDP no estacionaria 


El modelo más simple que se puede establecer 
para describir la tendencia de la serie anual de 
máximos es el logarítmico lineal con pendiente f 
estadísticamente diferente de cero. Este enfoque 
es adecuado en series de valores hidrológicos 
extremos, pues al emplear logaritmos de sus 
datos se puede aceptar que la muestra procede 


de una población normal. El modelo es (Vogel 
et al., 2011): 


y, =In(x,)=0+B-t+e, (2) 


siendo a. la ordenada al origen y e, los errores 
del modelo, los cuales proceden de las variacio- 
nes de las observaciones respecto a la recta de 
regresión; esta línea, por sí misma, proporciona 
la forma en que y, varía a través del tiempo, es 
decir: 


u, (1)=0+B: (3) 


Como la solución de mínimos cuadrados de 
los residuos para la ordenada al origen es (Ryan, 
1998): 


A 


a=y-p i (4) 


Entonces, sustituyendo la ecuación (4) en 


la (3) se obtiene el modelo de tendencia de la 
media: 


(0) =7+b- (1-7) 65) 


en donde f es la estimación de mínimos cuadra- 
dos de los residuos, cuya expresión es (Ryan, 
1998): 


p- La — (6) 


además, las medias son: 


_ 14 

y= 2 (7) 
- lu n+1 

A (8) 


Como ya se indicó, la ecuación (1) es un mo- 
delo probabilístico estacionario, pues asume que 
los momentos de y, dados por u, y G, no varían 
con el tiempo. Entonces, la sustitución del mo- 
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delo de tendencia de la media o ecuación (5) en 
la ecuación (1) conduce al modelo probabilístico no 
estacionario desarrollado y propuesto por Vogel 
et al. (2011): 


x,() exp [y+B-(1-5)+2,05, (9) 


en el cual 5, es la estimación de la desviación es- 
tándar de los logaritmos naturales de los datos 
hidrológicos extremos anuales y t es el contador 
de años, variando de 1 a n. En realidad, éste es 
el modelo más simple que se puede establecer 
para la tendencia, teniendo al tiempo como 
única covariable explicativa (López-De-La-Cruz 
éz Francés, 2014). S, se calcula con la expresión: 


(10) 


Factor de aumento de la predicción 


Como el propósito de los análisis es estimar el 
impacto en las predicciones de la tendencia ob- 
servada en el registro hidrológico de extremos 
anuales, primero se evaluará el denominado 
factor de aumento de la predicción (FAP), definido 
como el cociente entre la predicción estimada 
en un lapso futuro y la del presente; en relación 
con la ecuación (9), se obtiene como (Vogel et 
al., 2011): 


FAP = E =exp($-At) (11) 
x,(0) 

La expresión anterior indica que el incre- 
mento o decremento en una predicción en un 
tiempo Af del futuro sólo depende de la tenden- 
cia $ estimada en el registro, sea ésta positiva 
o negativa. Vogel et al. (2011) encuentran en 
cuencas localizadas en las áreas periféricas de 
las grandes ciudades de Estados Unidos, como 
son Los Ángeles, Chicago y Boston, que los 
FAP en una década a futuro y para registros de 
crecientes con tendencias ascendentes varían de 
2 a5; es decir, que detectaron pendientes severas 
fluctuando de 0.070 a 0.160. En cambio, para re- 
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gistros de crecientes con tendencia descendente, 
los factores de disminución de la predicción (FDP) 
de un decenio a futuro oscilaron de 0.94 a 0.68. 


Reducción del periodo de retorno 


La reducción de la recurrencia será el tiempo 
promedio entre dos predicciones: la del futuro 
con periodo de retorno T, ligado con el tiempo 
t,, y la actual con valores correspondientes 
T, y t, Cada predicción está vinculada con 
estas probabilidades de excedencia p,= 1 / T,y 
p,=1/T,, Igualando las predicciones, estimadas 
con la ecuación (9), se podrá obtener el 1, busca- 
do, esto es (Vogel et al., 2011): 


pieles, 5 


n+1 


pablo Les, 5, (12) 


=exp 


Despejando, se obtiene: 


en donde D(-) es función de densidad de proba- 
bilidad de una variable normal estandarizada y 
representa la probabilidad de no excedencia de 
que tal variable sea menor que el valor del pa- 
réntesis. En el análisis de registros con tendencia 
descendente, el T, aumentará. 


Prueba de la pendiente del modelo logarítmico 
lineal 


En relación con la ecuación (2), su establecimien- 
to implica procesar registros preferentemente sin 
datos anuales faltantes, pues es una regresión 
lineal; además, se debe probar que su pendiente 
B es estadísticamente diferente de cero (Kund- 
zewicz Robson, 2004; Machiwal €: Jha, 2008). 
Para ello se aplicó la prueba sugerida por Ostle 
y Mensing (1975) basada en el estadístico de la 
t de Student, que utiliza las varianzas de los 
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errores y de la pendiente de la recta de regresión 
(S? y S; ). Si el valor absoluto calculado para 
Et (ecuación (14)) es mayor que el crítico (Et), 
obtenido para la distribución t de Student con 
v =1 - 2 grados de libertad y a = 5%, en una 
prueba de dos colas, la pendiente Bes significa- 
tiva, esto es, existe tendencia lineal. El problema 
de tal prueba es que no distingue entre persis- 
tencia y tendencia (Adeloye €: Montaseri, 2002). 


Epa LB (14) 


siendo: 


S, e E== (15) 


A (16) 


En la expresión anterior, / es el valor esti- 
mado con la recta de regresión, es decir, con la 
ecuación (2). Para el cálculo del estadístico Et. 
se utilizó el algoritmo computacional propuesto 
por Zelen y Severo (1972), cuyas expresiones 
son: 


Et (17) 


=Z + gap) g(z,) ¿26 Ñ 842) 
v 


Py y? Y 
gr(2)=(2+2)/4 
gr (2) = (52 +16z% +32) /96 
ga(2)=(327 +192* +172? -152) / 384 
2,(2)= (792+ 7762 +14822*-1920z*- 9452) /92160 


siendo z, la variable normal estandarizada, 
cuyos valores para niveles de significancia (a) 
de 10, 5 y 1%, en una prueba de dos colas son 
1.64485, 1.95996 y 2.57583; v son los grados de 
libertad. 


Aproximaciones para la distribución 
normal 


En las ecuaciones (1) y (9) se requiere la eva- 
luación de la variable normal estandarizada 2 
la cual se estimó con la aproximación racional 
expuesta por Zelen y Severo (1972). Cuando la 
probabilidad de excedencia p varía de 0 a 0.50, 
sus expresiones son: 


3 
4,+4,'9+0,*g 


z,=9- 18 
e 1+b,:q+b,-q +b,:q* id 
siendo: 
q=yIn(1/p?) (19) 
y 
q, =2.515517 — a, =0.802853 a, = 0.010328 
b,=1.432788  b,=0.189269  b,=0.001308 


Cuando p > 0.50, se usa 1 - p en la ecuación 
(19) y se le cambia el signo al valor calculado de 
z, con la ecuación (18). 

Para estimar la probabilidad normal D(-) de 
no excedencia en la ecuación (13) se utilizó una 
aproximación polinomial también sugerida por 
Zelen y Severo (1972), cuya expresión es: 


D(arg) =1-Z(arg) 
(e, w+c, w+c,w+c,w+c,w) (20) 
en la cual: 
w=1/(1+d:arg) (1) 
d =0.2316419 


c, =-0.356563782 
c,=-1.821255978 


c, = 0.31938153 
e, =1.781477937 
c¿=1.330274429 


En la ecuación (20), Z(arg) es la ordenada de 
la curva normal calculada con la ecuación: 


li (22) 
TT 
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Error estándar de ajuste 


Desde mediados de la década de los años setenta 
se formuló el error estándar de ajuste (EEA) como 
una medida cuantitativa que estima la calidad 
estadística del modelo probabilístico ajustado 
y que además permite la comparación objetiva 
entre los diversos modelos que se prueban o 
ajustan a una muestra, al tener las unidades 
de los datos. Su expresión es la siguiente (Kite, 
1977): 


(23) 


en donde 1 es el número de datos de la muestra; 
x, los datos ordenados de menor a mayor; X, los 
valores estimados con la ecuación (9), en la cual 
t varía de lan y le va correspondiendo una pro- 
babilidad de no excedencia P(X < x), estimada 
con la fórmula de Weibull (Benson, 1962): 


MX<)- E (24) 
n+1 
en la cual m es el número de orden del dato, con 
1 para el menor y n para el mayor. 


Datos y resultados de las aplicaciones 
numéricas 


Estimaciones del análisis de frecuencias 


La descripción detallada del modelo probabilís- 
tico no estacionario de los incisos anteriores se 
puede resumir en las seis estimaciones siguien- 
tes: 1) se calculan los parámetros $ y á del mo- 
delo logarítmico lineal con las ecuaciones (6) y 
(4), y se prueba si el primero es estadísticamente 
diferente de cero por medio de las ecuaciones 
(14) a (17); 2) en seguida se calculan dos factores 
de aumento de las predicciones (FAP), uno para 
la primera década y otro para un lapso de 20 
años con base en la ecuación (11); 3) después se 
calculan dos reducciones del periodo de retorno 
de 100 años de una predicción a través de la 
ecuación (13), utilizando At de 10 y 20 años, 
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este cálculo y los dos siguientes utilizan las 
ecuaciones auxiliares (18) y (20); 4) haciendo 
uso de la ecuación (9), primero se calcula el 
error estándar de ajuste, según la ecuación (23); 
5) después se evalúan las curvas de cuantiles 
con periodos de retorno (Tr) de 2, 10, 50 y 100 
años, esto es, utilizando 2,= 0, 1.28155, 2.05375 
y 2.32635, respectivamente; 6) por último, se 
estiman las predicciones asociadas con los 
periodos de retorno de diseño empleados por 
lo común en los análisis de frecuencia de datos 
hidrológicos extremos, por medio de la ecuación 
(9) utilizando t = 1. 

Cuando la tendencia del registro procesado 
en ascendente, se estiman FAP y reducciones 
de la recurrencia de la predicción de Tr = 100 
años, así como predicciones mayores que las 
obtenidas con una FDP estacionaria. Cuando 
la tendencia es descendente, se estimarán FDP 
y aumentos en la recurrencia del evento de Tr 
= 100 años, además de predicciones menores 
que las obtenidas con una FDP estacionaria. 
Lo anterior se ilustra en las tres aplicaciones 
numéricas siguientes. 


Estación hidrométrica Badiraguato 
Crecientes anuales y su tendencia ascendente 


Esta estación de aforos pertenece a la Región 
Hidrológica núm. 10, Sinaloa, y se ubica en el río 
del mismo nombre. Su registro actualizado de 
gastos máximos anuales procede de la GASIR 
de la Comisión Nacional del Agua (Conagua), 
con clave 10079, abarca 40 años en el periodo de 
1960 a 1999, y se expone en la segunda columna 
del cuadro 1. Su análisis de homogeneidad de- 
tecta tendencia ascendente con las pruebas de 
Kendall y de Spearman (WMO, 1971; Kottegoda, 
1980; Kundzewicz € Robson, 2004; Machiwal éz 
Jha, 2008; Campos-Aranda, 2015). La aplicación 
de las ecuaciones (4) y (6) aporta Q = 5.8595 y 
$ = 0.023393. La pendiente estimada es impor- 
tante y, sin embargo, resulta estadísticamente 
igual a cero, ya que Et= 1.7697 < Et_= 2.0244; lo 
anterior debido a que la varianza de los errores 
(S?) resultó muy grande, en parte debido a la 
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enorme dispersión que muestran los datos al 
variar de 64 a 9 245 m*/s (ver columna 3 del 
cuadro 1). Al tomar en cuenta que el registro 
no presenta persistencia, pues su coeficiente 
de correlación serial de orden 1 (r,) es de sólo 
0.054, se decidió continuar con el ajuste del 
modelo probabilístico no estacionario; es decir, 
considerar que el registro de crecientes anuales 
de la estación de aforos Badiraguato presenta 
tendencia ascendente. 


Impacto de la tendencia ascendente en las 
predicciones 


Con base en la ecuación (11) y para dos lapsos 
a futuro, la primera y segunda décadas, se es- 
timó el factor de aumento de las predicciones 
en 1.264 y 1.597, lo cual significa que cualquier 
predicción estimada al final del registro en los 
siguientes 10 años aumentará un 26.4%, y en 
los subsiguientes 10 años llegará a casi un 60%. 
Por otra parte, y de acuerdo con la ecuación 
(13), la creciente de periodo de retorno de 100 
años estimada al final del registro Zo =2.32635) 
reducirá su recurrencia en la primera década a 
54.3 años y al término del segundo decenio a 
31.0 años. Las cantidades anteriores ponen de 
manifiesto la importancia del impacto de la ten- 
dencia ascendente observada en las crecientes 
de la estación Badiraguato en las predicciones. 


Error estándar de ajuste y curvas de cuantiles 


En la columna 4 del cuadro 1 se tienen las es- 
timaciones realizadas con la ecuación (9) para 
t, variando de l a n; z, asociada con la proba- 
bilidad definida con la fórmula de Weibull o 
ecuación (24), y S, = 0.978602. La aplicación de 
la ecuación (23) con los valores de las columnas 
3 y 4 del cuadro 1 conduce al error estándar de 
ajuste de 585 m*/s. De nuevo, la aplicación de la 
ecuación (9) permite estimar las curvas de cuan- 
tiles con periodos de retorno de 2, 10, 50 y 100 
años, cuyos valores se muestran en las columnas 
5 a 8 del cuadro 1. La curva de cuantiles de 
Tr = 2 años corresponde a la recta de regresión 
logarítmica, la cual se muestra en la figura 1 con 


los datos del registro de crecientes. Se observa 
en la figura 1 que arriba de la curva de cuantiles 
de Tr = 10 años quedan cuatro datos, lo cual es 
correcto, pues 1 = 40; también se detecta que las 
crecientes de los años 1962 y 1981 tienen un Tr 
superior a los 100 años. 


Predicciones con periodos de retorno de diseño 


La aplicación de la ecuación (9) con t=1 permite 
estimar las predicciones asociadas con los ocho 
periodos de retorno mostrados en el cuadro 2. 
Con objeto de verificar que las predicciones es- 
timadas con el modelo probabilístico no estacio- 
nario tienen magnitud mayor, se ajustaron a los 
datos las cuatro FDP más comunes en los análi- 
sis de frecuencias de crecientes, destacando que 
su aplicación no es válida, pues el registro de la 
estación Badiraguato mostró tendencia ascen- 
dente. Estas FDP se ajustaron con sus métodos 
más comunes, los cuales se pueden consultar 
en Stedinger, Vogel y Foufoula-Georgiou (1993); 
Hosking y Wallis (1997); Rao y Hamed (2000) 
y Campos-Aranda (2006). Las predicciones de 
las FDP aplicadas corresponden a su método de 
ajuste que condujo al menor EEA y se muestran 
en el cuadro 2. 

Se observa que el EEA obtenido con el mode- 
lo no estacionario es el menor y sus predicciones 
son superiores en los periodos de retorno bajos, 
y ello se cumple hasta el Tr de 100 años. Las FDP 
aplicadas conducen a predicciones muy grandes 
en los periodos de retorno elevados (Tr > 100 
años) debido a presencia de los valores extremos 
de los años 1981, 1996 y 1962 (ver cuadro 1). 


Estación hidrométrica Puente Subpacífico 
Crecientes anuales y su tendencia descendente 


Esta estación de aforos también pertenece a la 
Región Hidrológica núm. 10, Sinaloa. Mide las 
descargas del río Culiacán después de su paso 
por la ciudad del mismo nombre, capital del 
estado, y su área drenada es de 15 731 km?. Al 
río Culiacán lo forman dos grandes colectores: 
el río Humaya y el río Tamazula; en el segundo 
se localiza la presa Sanalona, construida durante 
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Cuadro 1. Gastos máximos anuales (m?*/s) en la estación hidrométrica Badiraguato y cálculos asociados con el error estándar 


de ajuste y con sus curvas de cuantiles. 


1 2 3 4 5 6 7) 8 
Datos Gasto estimado para: 
Año (t) 
Observado | Ordenado | Prob. Weibull | Tr=2 años Tr= 10 años Tr=50años | Tr=100 años 
1960 (1) 435 64 52.2 359 1 258 2 678 3 496 
1961 (2) 276 127 72.6 367 1 288 2741 3579 
1962 (3) 4220 213 90.7 376 1318 2 806 3 664 
1963 (4) 420 236 108.3 385 1349 2872 3751 
1964 (5) 236 242 126.0 394 1381 2 940 3 839 
1965 (6) 213 242 144.0 403 1 414 3010 3 930 
1966 (7) 982 248 162.8 413 1 447 3 081 4 023 
1967 (8) 480 250 182.4 428 1 482 3154 4 119 
1968 (9) 435 276 203.0 433 1517 3229 4216 
1969 (10) 64 277 224.7 443 1552 3305 4316 
1970 (11) 250 290 247.7 153 1589 3 383 4 418 
1971 (12) 435 305 27203 464 1627 3 464 4 522 
1972 (13) 470 316 298.4 475 1665 3 546 4 630 
1973 (14) 585 401 326.3 486 1705 3629 4 739 
1974 (15) 600 420 356.1 498 1745 3715 4851 
1975 (16) 290 435 388.2 510 1786 3803 4 966 
1976 (17) 316 435 422.7 522 1829 3 893 5 084 
1977 (18) 242 435 459.8 534 1872 3 985 5 204 
1978 (19) 248 461 499.9 547 1916 4 080 5 327 
1979 (20) 891 470 543.2 560 1962 4 176 5 453 
1980 (21) 461 480 590.3 573 2 008 4275 5 582 
1981 (22) 9 245 481 641.5 587 2 056 4.376 5714 
1982 (23) 2 388 487 697.4 600 2104 4 480 5 850 
1983 (24) 27) 488 758.7 615 2154 4 586 5 988 
E 1984 (25) 481 582 826.0 629 2205 4 695 6 130 
5 1985 (26) 614 585 900.4 644 22 4 806 6275 
E 1986 (27) 1064 600 982.8 659 231 4919 6 423 
S 1987 (28) 401 614 1 074.7 675 2.365 5 036 6576 
3 1988 (29) 305 891 1177.8 691 2421 5 155 6731 
El 1989 (30) 582 904 1 294.3 707 2 479 5277 6890 
a 1990 (31) 2 035 933 1 427.1 724 2 537 5 402 7054 
le 1991 (32) 488 982 1 580.0 741 2 597 5 530 7221 
E 1992 (33) 242 1 064 1 758.4 759 2 659 5661 7391 
5 1993 (34) 487 1193 1 969.9 777 2722 DOS 7566 
S 1994 (35) 933 2 035 2 226.2 795 2786 5 932 7745 
al 1995 (36) 127 2 388 2 545.6 814 2 852 6 072 7929 
S 1996 (37) 4 490 2529 2 960.5 833 2920 6216 8 116 
3 1997 (38) 904 4220 3 534.2 853 2 989 6 363 8 309 
8 1998 (39) 2529 4 490 4 418.5 873 3059 6514 8 505 
g 1999 (40) 1193 9 245 6 148.8 894 3132 6671 8710 
E 
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Figura 1. Diagrama de datos y cuantiles estimados con el modelo logarítmico lineal en la estación hidrométrica Badiraguato de 


la Región Hidrológica núm. 10, Sinaloa. 


Cuadro 2. Predicciones de diseño (m?*/s) en la estación hidrométrica Badiraguato estimadas con el modelo probabilístico no 


estacionario y su contraste con las obtenidas ajustando cuatro FDP comunes en el análisis de frecuencias de crecientes. 


Modelo probabilístico ajustado y EEA Periodos de retorno de diseño en años 

(método de ajuste) (m/s) 2 5 10 25 50 100 500 1 000 
Ny Ea nonanO pepresión esa 585 | 894 | 2036 | 3132 | 4958 | 6671 | 8710 | 14946 | 18392 
lineal) 
e tipo III (momentos en dominio [2 [487 | 1200 | 2122 | 4215 | 6870 | 11.004 | 31415 | 48707 
General de valores extremos (sextiles) 765 499 1234 | 3067 | 3799 | 5864 8944 | 23300 | 35 016 
Pareto generalizada (momentos L) 808 472 1143 | 1946 | 3639 | 5665 8683 | 22674 | 34 027 
Logística generalizada (momentos L) 888 499 1104 | 1812 | 3353 | 5282 8293 | 23564 | 36 946 


los años 1940 a 1948, con un área de cuenca de 
3 250 km?. En el río Humaya se ubica la presa 
Adolfo López Mateos, construida durante los 
años 1957 a 1964, con un área de cuenca de 
10 972 km? (Aldama, Ramírez, Aparicio, Mejía- 
Zermeño, £ Ortega-Gil, 2006). De acuerdo con 
los datos anteriores, la estación hidrométrica 
Puente Subpacífico tiene una cuenca libre de 
sólo 1 509 km?, por lo cual su régimen hidro- 
lógico estará sumamente afectado por los dos 
aprovechamientos descritos. 


Su registro de crecientes anuales procede 
del sistema BANDAS (IMTA, 2003) y abarca del 
año 1924 a 1992, contando con 66 datos, pues 
el lapso de 1959 a 1961 está faltante. De modo 
sorpresivo, la creciente máxima registrada de 
11 000 m?/s ocurrió en 1943, durante el periodo 
de construcción de la presa Sanalona, y la última 
creciente importante aconteció en 1958, con 
5 715 m?/s, ya iniciado el proceso constructivo 
de la presa Adolfo López Mateos, cuyo impacto 
cambia de forma notable el régimen hidrológico, 
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como se observa en la figura 2, a partir del dato 
número 36. Para este registro, en el cuadro 3, 
similar al cuadro 1, se citan en exclusiva sus va- 
lores y estimaciones por quinquenios. El análisis 
de homogeneidad del registro citado detecta 
persistencia (r, = 0.366), y tendencia descenden- 
te con las pruebas de Kendall y de Spearman 
(WMO, 1971; Kottegoda, 1980; Kundzewicz éz 
Robson, 2004; Machiwal é Jha, 2008; Campos- 
Aranda, 2015). La aplicación de las ecuaciones 
(4) y (6) aporta ú = 8.1150 y $ = -0.027661. 
La pendiente estimada es importante y por ello 
resulta estadísticamente diferente de cero, ya 
que Et=|-6.4115| > Et, = 1.9977. Un valor de 
0.842596 resultó para S,, 


Impacto de la tendencia en las predicciones 


Con base en la ecuación (11) se estimó el factor 
de disminución de las predicciones (FDP) en 
0.758 y en 0.575, para la primera y segunda 
décadas a futuro; son reducciones importantes 
debido a que la pendiente f negativa estimada 
es fuerte. En relación con la creciente de periodo 
de retorno 100 años calculada al final del regis- 


tro, de acuerdo con la ecuación (13), aumentará 
su recurrencia en la primera década a 252.2 años 
y en la segunda década a 701.4 años. 


Error estándar de ajuste y curvas de cuantiles 


Siguiendo un procedimiento idéntico al descrito 
para este mismo inciso, se obtuvo un error es- 
tándar de ajuste de 1 652 m* /'s, así como las cur- 
vas de cuantiles con periodos de retorno de 2, 
10, 50 y 100 años, cuyos valores por quinquenios 
se muestran en las columnas 5 a 8 del cuadro 
3 y en la figura 2, con los datos del registro de 
crecientes. La curva de cuantiles de Tr = 2 años 
corresponde a la recta de regresión logarítmica. 


Predicciones con periodos de retorno de diseño 


La aplicación de la ecuación (9) con t = n per- 
mite estimar las predicciones asociadas con 
los ocho periodos de retorno mostrados en el 
cuadro 4. La verificación de la magnitud menor 
de estas predicciones se realizó aplicando las 
cuatro FDP más utilizadas en los análisis de 
frecuencias de crecientes, ajustada cada una 


Cuadro 3. Gastos máximos anuales (m?/s) en la estación hidrométrica Puente Subpacífico y cálculos asociados con el error 


estándar de ajuste y con sus curvas de cuantiles expuestos por quinquenios. 


1 2 3 4 5 6 7 8 
O Datos Gasto estimado para 
Observado | Ordenado | Prob. Weibull | Tr=2años | Tr=10años | Tr=50años | Tr=100 años 
1924 (1) 4.562 330 521.6 3253 9 578 18 358 23 099 
1928 (5) 1122 499 863.8 2 MZ 8574 16 435 20 679 
1933 (10) 1780 612 1 056.2 2 536 7 467 14 313 18 008 
1938 (15) 2412 651 1 165.2 2 209 6502 12 464 15 682 
1943 (20) 11 000 710 1 232.3 1923 5 663 10 854 13 657 


1948 (25) 790 852 1275 1675 4931 9 452 11 893 
1953 (30) 3 490 944 1 306.1 1 459 4 294 8 231 10 357 
1958 (35) Dl 77) 1331.5 1270 3740 7 168 9019 
1966 (40) 944 1173 1 360.0 1106 3257 6242 7854 
1971 (45) 900 2335 1 400.3 963 2 836 5 436 6 839 
1976 (50) 1090 2760 1 465.9 839 2470 4 734 5956 
1981 (55) 1173 9/99 1584.1 731 2151 4 122 5 187 
1986 (60) 685 4755 1 833.9 636 1873 3590 4 517 
1992 (66) 651 11 000 3 360.9 539 1587 3 042 3 827 
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Figura 2. Diagrama de datos y cuantiles estimados con el modelo logarítmico lineal en la estación hidrométrica Puente 
Subpacífico de la Región Hidrológica núm. 10, Sinaloa. 


Cuadro 4. Predicciones de diseño (m?/s) en la estación hidrométrica Puente Subpacífico estimadas con el modelo probabilístico 
no estacionario y su contraste con las obtenidas ajustando cuatro FDP comunes en el análisis de frecuencias de crecientes. 


Modelo probabilístico ajustado y EEA Periodos de retorno de diseño en años 
(método de ajuste): (m/s) 2 5 10 25 50 100 500 1000 


No estacionario (regresión logarítmica 


. 1 652 539 1095 | 1587 | 2356 | 3042 | 3827 | 6093 | 7284 
lineal) 


Pareto generalizada (momentos L) 461 1253 | 2777 | 4154 | 6330 | 8296 | 10588 | 17486 | 21 302 


Log-Pearson tipo III (momentos en dominio [— ¿,, |, 369 | 2915 | 4198 | 6047 | 7558 | 9159 | 13198 | 15064 


real) 

General de valores extremos (momentos L) 516 1319 | 2643 | 3893 | 6081 8300 | 11177 | 21 645 | 28 545 
Logística generalizada (momentos L) 549 1331 | 2593 | 3777 | 5905 | 8145 | 11159 | 22898 | 31 130 
con sus métodos más comunes (Stedinger et al., FDP empleadas son bastante flexibles y se 
1993; Hosking € Wallis, 1997; Rao €: Hamed, ajustan mejor a los datos dispersos o extremos 
2000, Campos-Aranda, 2006); los resultados se del registro; pero se aclara que su aplicación no 
exponen en el cuadro 4 para los EEA menores. es válida, pues el registro de crecientes anuales 
Se observa que las FDP reportan EEA mucho de la estación de aforos Puente Subpacífico es 

más bajos, debido principalmente a que las no estacionario. 


Tecnología y Ciencias del Agua, vol. VII, núm. 3, mayo-junio de 2016, pp. 171-186 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VII, núm. 3, mayo-junio de 2016, pp. 171-186 


Campos-Aranda, Modelo probabilístico simple para análisis de frecuencias en registros hidrológicos extremos con tendencia 


Estación pluviométrica Zacatecas 
PMD anual y su tendencia ascendente 


El registro de precipitación máxima diaria 
(PMD) anual de la estación climatológica 
Zacatecas fue proporcionado y verificado en 
sus valores extremos por la Dirección Local 
Zacatecas de la Conagua y consta de 58 datos 
en el lapso de 1953 a 2010; sus valores se expo- 
nen en la columna 2 del cuadro 5. Su estudio de 
homogeneidad detecta tendencia ascendente 
con las pruebas de Kendall y de Spearman 
(WMO, 1971; Kottegoda, 1980; Kundzewi- 
cz € Robson, 2004; Machiwal éz Jha, 2008; 
Campos-Aranda, 2015) y no tiene persistencia 
(r, = 0.002). Los parámetros del modelo loga- 
rítmico lineal, según las ecuaciones (6) y (4) 
resultaron ser á = 3.6287 y f$ = 0.005679. Con 
base en las ecuaciones (14) a (17), se obtuvo un 
Et = 2.3361 > Et = 2.0032, de manera que esta 


prueba indica que la pendiente ascendente B 
es estadísticamente diferente de cero. Además, 
S, = 0.319011. 


Impacto de la tendencia en las predicciones 


Con base en la ecuación (11) se estimó el factor 
de aumento de las predicciones (FAP) en un 
5.8% y en un 12%, para la primera y segunda 
décadas a futuro; son incrementos reducidos de- 
bido a que la pendiente $ estimada es realmente 
baja. En relación con la creciente de periodo de 
retorno 100 años estimada al final del registro, 
de acuerdo con la ecuación (13) reducirá su 
recurrencia en la primera década a 63.2 años y 
en la segunda década a 41.0 años. 


Error estándar de ajuste y curvas de cuantiles 


A partir de un procedimiento idéntico al descri- 
to para este mismo inciso, se obtuvo un error es- 


Cuadro 5. PMD anual (mm) en la estación pluviométrica Zacatecas y cálculos asociados con el error estándar de ajuste 


y sus curvas de cuantiles. 


1 2 3 4 5 6 Y 8 
Datos PMD estimada para: 
Año (t) 
Observado | Ordenado | Prob. Weibull | Tr=2 años Tr= 10 años Tr=50 años | Tr=100 años 
1953 (1) 35.1 16.4 19.3 37.9 57.0 72.9 79.6 
1954 (2) 23.0 23.0 21.3 38.1 57.3 73.3 80.0 
1955 (3) 82.5 26.4 PL 38.3 57.7 73.8 80.5 
1956 (4) 64.0 27.0 23.9 38.5 58.0 74.2 80.9 
1957 (5) 36.9 28.9 25.0 38.7 58.3 74.6 81.4 
1958 (6) 43.5 29.0 26.0 39.0 58.7 75.0 81.9 
1959 (7) DRA po 26.9 39.2 59.0 SS 82.3 
1960 (8) 38.2 30.0 27.7 39.4 59.3 75.9 82.8 
1961 (9) 30.0 32.0 28.6 39.6 59.7 76.3 83.3 
1962 (10) 26.4 32.0 29.4 39.9 60.0 76.8 83.7 
1963 (11) 28.9 35.0 30.2 40.1 60.3 pz 84.2 
1964 (12) 38.0 35.1 30.9 40.3 60.7 77.6 84.7 
1965 (13) 39.4 35.4 31.7 40.5 61.0 78.1 85.2 
1966 (14) 39.2 36.0 32.5 40.8 61.4 78.5 85.7 
1967 (15) 46.4 36.9 33.2 41.0 61.7 79.0 86.1 
1968 (16) 48.4 38.0 34.0 41.2 62.1 79.4 86.6 
1969 (17) 16.4 38.2 34.7 41.5 62.4 79.9 87.1 
1970 (18) 44.4 39.2 35.5 41.7 62.8 80.3 87.6 
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Cuadro 5 (continuación). PMD anual (mm) en la estación pluviométrica Zacatecas y cálculos asociados con el error estándar 


de ajuste y sus curvas de cuantiles. 


1 2 3 4 5 6 7 8 
AñO(D Datos PMD estimada para: 
Observado | Ordenado | Prob. Weibull | Tr=2 años Tr= 10 años Tr=50 años | Tr=100 años 
1971 (19) 48.4 39.4 36.2 42.0 63.1 80.8 88.1 
1972 (20) 44.0 40.0 37.0 42.2 63.5 81.2 88.6 
1973 (21) 75.0 40.0 37.7 42.4 63.9 81.7 89.1 
1974 (22) 29.0 40.9 38.5 42.7 64.2 82.2 89.6 
1975 (23) 57.0 41.0 39.3 42.9 64.6 82.6 90.1 
1976 (24) 36.0 41.2 40.0 43.2 65.0 83.1 90.7 
1977 (25) 50.0 43.5 40.8 43.4 65.3 83.6 2 
1978 (26) 35.4 43.5 41.6 43.7 65.7 84.1 91.7 
1979 (27) 40.0 44.0 42.4 43.9 66.1 84.5 9212 
1980 (28) 32.0 44.0 43.3 44.2 66.5 85.0 92.7 
1981 (29) TD 44.4 44.1 44.4 66.8 85.5 93.3 
1982 (30) 29.2 44.5 45.0 44.7 67.2 86.0 93.8 
1983 (31) 43.5 46.0 45.8 44.9 67.6 86.5 94.3 
1984 (32) 50.0 46.4 46.7 45.2 68.0 87.0 94.9 
1985 (33) 60.3 46.5 47.6 45.4 68.4 87.5 95.4 
1986 (34) 75.0 48.4 48.6 45.7 68.8 88.0 96.0 
1987 (35) 52.5 48.4 49.5 45.9 69.1 88.5 96.5 
1988 (36) 46.0 50.0 50.5 46.2 69.5 89.0 97.1 
1989 (37) 40.0 50.0 51.5 46.5 69.9 89.5 97.6 
1990 (38) 65.7 51.0 52.6 46.7 70.3 90.0 98.2 
1991 (39) 52.0 SL 53.7 47.0 70.7 90.5 98.7 
1992 (40) 32.0 52.0 54.8 47.3 71.1 91.0 99.3 
1993 (41) 44.5 52.0 55.9 47.5 ZALAS 91.5 99.8 
1994 (42) 54.0 52.2 57.1 47.8 72.0 92.1 100.4 
1995 (43) 46.5 525 58.4 48.1 72.4 92.6 101.0 
1996 (44) 63.0 53.0 59.7 48.4 72.8 93.1 101.6 
1997 (45) 40.9 54.0 61.1 48.6 Sa 93.6 102.1 
1998 (46) 67.5 57.0 62.5 48.9 73.6 94.2 102.7 
1999 (47) 41.2 60.3 64.1 49.2 74.0 94.7 103.3 
2000 (48) 35.0 63.0 65.7 49.5 74.4 95.2 103.9 
2001 (49) 51.5 64.0 67.5 49.7 74.9 95.8 104.5 
2002 (50) 41.0 64.0 69.4 50.0 75.3 96.3 105.1 
2003 (51) 44.0 65.0 ZALOS 50.3 DI 96.9 105.7 
2004 (52) 52.0 65.7 73.8 50.6 76.2 97.4 106.3 
2005 (53) 27.0 67.5 76.4 50.9 76.6 98.0 106.9 
2006 (54) 53.0 71.0 79.3 51.2 77.0 98.5 107.5 
2007 (55) 64.0 TES 82.9 LS RS 99.1 108.1 
2008 (56) 71.0 75.0 87.3 51.8 77.9 99.7 108.7 
2009 (57) 51.0 75.0 93.2 52.1 78.4 100.2 109.3 
2010 (58) 65.0 82.5 103.0 52.4 78.8 100.8 110.0 
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Figura 3. Diagrama de datos y cuantiles estimados con el modelo logarítmico lineal en la estación hidrométrica Zacatecas de la 


Región Hidrológica núm. 37, El Salado. 


Cuadro 6. Predicciones de diseño (mm) en la estación pluviométrica Zacatecas estimadas con el modelo probabilístico no 


estacionario y su contraste con las obtenidas, ajustando cuatro FDP comunes en el análisis de frecuencias de crecientes. 


Modelo probabilístico ajustado y EEA Periodos de retorno de diseño en años 
(método de ajuste) (mm | 2 5 10 | 25 | 50 | 100 | 500 | 1000 
No estacionario (regresión logarítmica lineal) 5.7 52 68 79 92 101 110 | 131 140 
General de valores extremos (sextiles) 15 45 58 66 US 81 87 | 98 102 
Log-pearson tipo IT (momentos en dominio real) 17 45 58 66 75 81 87 | 99 104 
Logística generalizada (momentos L) 1.8 45 57 65 75 83 92) | 114 125 
Pareto generalizada (momentos) 2.3 45 61 68 74 77 79 [os 81 


tándar de ajuste de 5.7 milímetros, así como las 
curvas de cuantiles con periodos de retorno de 
2, 10, 50 y 100 años, cuyos valores se muestran 
en las columnas 5 a 8 del cuadro 3, y en la figura 
3 con los datos del registro. Se insiste en que la 
curva de cuantiles de Tr =2 años corresponde a 
la recta de regresión logarítmica. Se observa en 
la figura 3 que arriba de la curva de cuantiles 
de Tr = 10 años quedan cinco datos, lo cual es 
correcto, pues n = 58; también detecta que la 
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PMD máxima del año 1955 tiene un Tr del orden 
de 100 años. 


Predicciones con periodos de retorno de diseño 


Por último, la aplicación de la ecuación (9) con 
t = 1 permite estimar las predicciones asociadas 
con los ocho periodos de retorno mostrados en 
el cuadro 6. La verificación de la magnitud mayor 
de estas predicciones se realizó aplicando las 
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cuatro FDP más utilizadas en los análisis de 
frecuencias de crecientes, ajustada cada una 
con sus métodos más comunes (Stedinger et al., 
1993; Hosking € Wallis, 1997; Rao €z Hamed, 
2000, Campos-Aranda, 2006); los resultados se 
exponen en el cuadro 6 para los EEA menores. 
Se observa que las FDP reportan EEA más bajos, 
debido en parte a que este registro no presenta 
valores dispersos e incluso las FDP ajustadas 
muestran límite superior; pero se insiste que su 
aplicación no es válida, pues el registro de PMD 
anual de la estación pluviométrica Zacatecas es 
no estacionario. 


Conclusiones 


El modelo probabilístico no estacionario desa- 
rrollado y propuesto por Vogel et al. (2011), para 
los análisis de frecuencias, descrito con detalle 
en este trabajo, es quizás el enfoque más simple 
y práctico que toma en cuenta la tendencia 
(ascendente o descendente) encontrada en 
ciertos registros de datos hidrológicos extremos, 
pues considera como única covariable explicati- 
va al tiempo, con respecto al cual varía la media 
del modelo ajustado, que es la distribución log- 
normal de dos parámetros. 

Dos aspectos prácticos importantes del 
modelo son las estimaciones de los factores de 
aumento de las predicciones (ecuación 11) y 
de la reducción de la recurrencia de una cierta 
predicción (ecuación 13), pues permiten estimar 
numéricamente el impacto real de la tendencia 
ascendente observada en las predicciones, como 
al final de una década o al término de 20 años 
hacia el futuro. 

Cuando un registro de datos hidrológicos ex- 
tremos muestre tendencia (ascendente o descen- 
dente) y ésta resulte estadísticamente diferente 
de cero, un primer modelo probabilístico que 
se recomienda probar para abordar el análisis 
de frecuencias no estacionario es el descrito por la 
ecuación (9), en particular debido a su extraor- 
dinaria sencillez de uso e interpretación de sus 
diferentes resultados, como se mostró en las tres 
aplicaciones numéricas descritas con detalle. 
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